Soliton Instabilities and Vortex Streets Formation in a Polariton Quantum Fluid 
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Exciton-polaritons have been shown to be an optimal system in order to investigate the properties 
of bosonic quantum fluids. We report here on the observation of dark solitons in the wake of 
engineered circular obstacles and their decay into streets of quantized vortices. Our experiments 
provide a time-resolved access to the polariton phase and density, which allows for a quantitative 
study of instabilities of freely evolving polaritons. The decay of solitons is quantified and identified 
as an effect of disorder-induced transverse perturbations in the dissipative polariton gas. 



Perturbations of quantum fluids can lead to the cre- 
ation of solitary waves called solitons resulting from the 
compensation between dispersion and particle interac- 
tion. In the particular case of repulsive interaction, dark 
solitons are created. These density depressions move in 
the fluid while keeping a constant shape and they are 
characterized by a phase jump across the density mini- 
mum. Since the first theoretical prediction [l|, dark soli- 
tons have been studied and then observed in a variety 
of systems such as nonlinear lattices [2], thin magnetic 
films [3] and complex plasma [4]. They have attracted 
considerable interest especially in the field of nonlinear 
optics [H[ , because of their consequent use in communica- 
tion devices (i.e optical fibers Q), and in atomic BEC Q- 
As quantized vortices (8l4lo| , dark solitons are BEC exci- 
tations, which arise spontaneously upon the phase tran- 
sition. As such, they are clear evidences for the onset of 
a quantum behavior and powerful tools to understand 
BEC instabilities. Controlling these instabilities is of 
crucial importance for the development of optoelectronic 
devices based on quantum fluids in which stable regimes 
and structures are required. Dark solitons are considered 
as the dispersive and nonlinear analog of shock waves of 
supersonic motion 11]. The creation of solitons by phase 
imprinting in BEC has been reported 12, 13j |. triggering 



a growing interest in their hydrodynamic formation and 
stability of solitons. 

In this paper we report on the observation of hy- 
drodynamic oblique dark solitons in a 2D polariton 
fluid and the formation of quantized vortex streets. 
Polaritons are bosonic quasi-particles arising from the 
strong coupling between quantum well excitons and 
photons in semiconductor microcavities. Due to their 
small effective mass and to their strong nonlinearity, 
they have turned out to be an optimal system in order 
to investigate the properties of a bosonic quantum fluid. 
Polaritons can undergo BEC [14[ and, by virtue of their 
non-equilibrium nature, they are accompanied with the 
spontaneous appearance of quantized vortices ll5| . With 



the hydrodynamic properties of a polariton superfluid 
have been investigated both theoretically 18|, [l9j and 
experimentally. Nucleation of quantized vortex pairs [lo[ 
and dark solitons 21] have been observed in the wake of 
an obstacle inserted into the polariton fluid. Differently 
to the work reported in Ref. [2l| , in which a steady state 
was created by a continuous wave pumping, our pulsed 
experiments allow us to retrieve the dynamics of the 
polariton fluid and thus to reconstruct the temporal 
behavior of dark solitons. Time and phase resolution 
consents to observe the formation of vortex streets 
arising from the unstable nature of solitons and to 
quantitatively access the stability conditions. 
It is well known that dark solitons are unstable with 
respect to transverse perturbations 22|, Q and that 



the demonstration of polariton superfluidity [16|, [17] 
much effort has been performed to better understand 
the nature of different turbulences arising from the 
breakdown of this fascinating state of matter. Recently, 



they eventually decay into other more stable structures 
[UQ- The driven dissipative nature of polaritons 
results in the exploration of different hydrodynamic 
regimes enriching the achievable features of polariton- 
solitons. This allows in particular the study of soliton 
instability. 



In our experiments we observe solitons in the wake 
of an engineered potential. The problem of a quantum 
fluid passing an obstacle can be formulated using the 
Gross-Pitaevskii equation and introducing an external 
potential term [8j. In the case of a polariton gas, whose 
nature is strongly dissipative due to the finite lifetime, 
it has recently been proved that the ratio between 
the local flow velocity and local value of the speed of 
sound, namely the Mach number M, determines the 
hydrodynamic regimes [l8[. Solitons are predicted to 
appear for M larger than 1 with values depending on the 
nature of the obstacle. We resonantly inject polaritons 
in a GaAs microcavity with a single InGaAs quantum 
well featuring a 3.5 meV Rabi splitting [27j. The initial 
in-plane momentum is imposed by the pump. A 2 ps 
long pulse is used to create a polariton wavepacket which 
can subsequently evolve freely in the microcavity plane 
during the polariton lifetime and then show features 
typical of quantum fluid hydrodynamics 120] . In the 
microcavity a set of different engineered circular traps, 
called mesas, have been realized by etching the spacer 
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FIG. 1: Time evolution of density (upper panel) and phase (lower panel) of a polariton wavepacket scattering on an engineered 
circular obstacle (green circles) of 3 \im diameter. Polaritons are injected at t — ps with a 1.5 \irn~ 1 initial momentum and a 
pump power of 5 mW in the vicinity of the mesa (column I). Solitons (II) are visible few picoseconds later and are characterized 
by low density straight lines and a phase shift. The soliton decay (III) comes along with the breaking of the phase- fronts and 
the formation of vortex streets (blue circles). The latter is pointed out by the apparition of several density minima coinciding 
with phase singularities. The motion of vortices is tracked along the flow during polariton lifetime (IV). 



before growing the top DBR [27[ . They result in negative 
potentials barriers for the resonantly injected polariton 
wavepacket. Phase and amplitude of polaritons are 
resolved in time through homodyne detection. This is 
achieved by analyzing the spatial interference between 
a local oscillator and the polariton emission out of a 
Mach-Zehnder interferometer [20|, [28|. The dynamics is 
reconstructed from the subsequent interference images 
taken by varying the delay between the two arms of the 
interferometer. 



Figure 1 shows an example of the scattering dynamics 
of a supersonic polariton wavepacket against a 3 \im 
mesa (green circles). Polaritons are injected with a 1.5 
//m -1 initial momentum and a pump power of 5 mW 
distributed over a laser spot of about 20 jim diameter 
(in all panels flow is from the bottom to the top). These 
parameters unambiguously place the experiment in the 
supersonic regime [20[. The FWHM and the position of 
the laser spot with respect to the mesa is indicated by 
a white dashed line in the density profile of Fig. 1-1. It 
can be observed in Fig.l-II that, due to the attractive 
nature of the obstacle, the injection pulse also excites 
confined polariton states close to the continuum 28, 29| . 
which are characterized by bright lobes. 
Hydrodynamic soliton formation is the result of the 
interaction between the polariton fluid and the mesa, 
which provides a continuous perturbation on the con- 
densate in motion. The dynamics of a 2D fluid passing 
an obstacle and soliton formation can be analyzed in 



both lab and fluid reference frame. In the latter, the 
mesa moves with respect to the fluid leaving behind 
perturbations which expand as circular waves. Positive 
interference between them occurs tangentially resulting 
in the oblique soliton formation. They then move with 
respect to the fluid frame in the y s direction whereas 
they grow in x s direction in the lab frame (see Figl-I). 
Such solitons are emphasized in Fig. l-II by green 
dashed lines. Phase jumps across the amplitude mini- 
mum can be observed in the phase map. 
Few picoseconds after the polariton injection, soli- 
tons decay into vortex streets. The breaking of the 
phase-fronts creates singularities associated with density 
minima: a clear evidence for quantized vortices. These 
are highlighted in the phase-map by blue and red circles. 
Column III shows the density profile and the phase map 
during the formation of vortex-antivortex pairs along 
the streets. The motion of vortices can be tracked in 
order to retrieve their dynamics (Fig. l-III and 1-IV). 



A unique way to investigate the soliton decay into vor- 
tex streets is given by the observation of the density 
changes in time along the soliton direction (x s in Fig. 
1-1). A plot versus time of the density map along x s is 
shown in Figure 2. After the pulse arrival at t = ps, 
solitons nucleate until t = 4 ps and their presence is 
revealed by a long low density region. The insert of Fig- 
ure 2 shows the characteristic soliton profile (black line) 
along its transversal direction y s at t = 3 ps, before the 
formation of vortices. During the whole soliton lifetime, 
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FIG. 2: Density map along the soliton direction x s versus 
time. After the pulse arrival at t — ps, the presence of a low 
density region stresses the nucleation of a soliton. Around 
t = 5 ps soliton decay occurs and the formation of vortices 
is highlighted by the appearance of periodic modulation of 
the density in the region previously occupied by the soliton. 
Insert plot shows the density (black) and the phase (blue) 
profile along y s at t = 3 ps calculated at around 7 \im from 
the mesa center showing the characteristic soliton transversal 
shape and the phase jump across the minimum. 



the soliton width has an almost constant value around 
4.5 /im, not limited by the healing length and indepen- 
dent of the laser pump power. The phase jump along y s 
(namely across the soliton lines) typical of dark solitons 
is shown in the insert (blue line). 

Around five picoseconds after the nucleation, solitons de- 
cay into vortex streets which are visualized by the ap- 
pearance of periodic modulations of the density. The 
presence of vortices in Fig. 2 is revealed by sloping stripes 
of alternate density. Low density valleys represent the 
motion of the cores of quantized votices. The size of the 
latter is estimated to be around 2 jim. Once formed, 
vortices move with a constant velocity v = 0.85 jam/ps 
along the Xg axis, as demonstrated by the slope of the 
core trajectories in Fig. 2. This value must be compared 
with the projection of flow velocity over x s . The latter is 
calculated considering the ballistic propagation of polari- 
tons with high momentum of injection [30| and it matches 
with the one extracted from the phase map at the time of 
injection. The aperture angle of an oblique soliton with 
respect to the flow direction is measured to be ~ 30° at 
the decay time giving a component of Vfi ow along x s of 
0.95 jim/ps. As expected, vortices move along the flow 
direction with a velocity comparable to that of the fluid. 

The same experiment is repeated for a set of excitation 
powers ranging from 0.5mW to 5mW. Contrarily to what 
could come out from a naive picture, we observe (see Fig. 
3a) that the solitons lines become less stable when the 
excitation power is increased (snapshots of the dynamics 
lower pump powers can be seen in the supplementary ma- 
terials). The black dots in the Fig. 3a show the average 



time delay for the vortex street formation. The hori- 
zontal bars represent the time window during which the 
soliton instability is observed, namely the time between 
the formation of the first and the last of the vortex pairs. 
Note that we consider the whole instability with vortices 
appearing progressively all along the solitons in region 
where the density and thus the fluid behavior could be 
at different stage of the soliton evolution. 
Our observations are in agreement with previous studies 
on hydrodynamic nucleation of solitons and can be un- 
derstood in terms of values of the Mach number in the ob- 
stacle perimeter. Solitons become more stable for higher 
values of M [Tl|, [HI. These are obtained experimentally 
by decreasing the sound speed, namely decreasing the 
polariton fluid density. 

Studies on ID solitons in atomic condensates revealed 
that the soliton speed v s is related to the ratio between 
the soliton amplitude and the fluid density n s /no through 
the relation v s /c s = yl — n s /no (no and n s are mea- 
sured as illustrated in the insert of Fig. 2 and c s is the 
speed of sound) [l3[. This means that v s basically van- 
ishes when the density at the bottom of the soliton line 
becomes very small. 

The time evolution of n s /no is shown in Fig. 3b. It 
has been calculated by integrating the polariton emis- 
sion over a 4 fim long area at around 7 fim away from 
the mesa. The Mach number is estimated from the val- 
ues of n s /no to be larger than 2, which confirms that our 
experiments are placed in the supersonic regime. During 
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FIG. 3: a - Pump power versus soliton lifetime. Horizontal 
bars show the time window at which vortex street formation 
occurs. Central black spots highlight that soliton lifetime in- 
creases while decreasing the pump power, b - Time evolution 
of the soliton amplitude n s over polariton fluid density no for 
different pump powers. Maximum of the ratio is found in cor- 
respondence of the soliton decay into vortex streets (dashed 
vertical lines). 




FIG. 4: Snapshots of the numerical simulations, based on 
the generalized Gross-Pitaevskii equations for the nucleation 
and the decay of dark solitons. The 2 fim diameter negative 
potential of —4 meV is indicated by a green cirlce. Flow 
direction is upstream. Random potential disorder is added in 
order to perturb solitons. 



the first picoseconds of the evolution, n s /no increases for 
all pump powers, showing that the solitons are becom- 
ing deeper. Remarkably, the observed break-up of the 
solitons corresponds quite precisely to the moment when 
this ratio reaches a critical value of 0.9 for all experimen- 
tal conditions. In disorder-free media, soliton instability 
is predicted to appear when n s /no reaches 1. It has also 
been shown that potential fluctuations modulate the fluid 
density and thus the local soliton velocity 32 , |33| . Fi- 
nally, in the case of soliton-like bubbles [34[, the critical 
ratio for stability is found to be smaller than 1 , depend- 



ing on the surrounding potential. From these theoretical 
inputs, we conclude that the critical ratio of 0.9 in our 
experiments can be attributed to the photonic disorder, 
typical of semiconductor microcavities. 

Our observations are supported by theoretical investi- 
gations performed by solving iteratively the generalized 
dissipative Gross-Pitaevskii equation for the lower polari- 
ton branch [35[ . Details about the methods can be found 
in Ref.(20[. Solitons nucleation in the wake of a negative 
potential and the subsequent decay into vortex streets is 
very well reproduced by the numerical simulations which 
take into account the pulsed excitation, the finite spot 
size of the laser, the exponential decay of the polariton 
population and the disorder potential. Snapshots of the 
simulated dynamics are reported in Fig. 4. They repro- 
duce perfectly the timescale of the experiments. After 
the polariton injection at t = ps, solitons are created al- 
most instantly and transform in vortex streets after about 
t = 5 ps. After creation, vortices move according to the 
flow as observed in the experiments. It is interesting to 
notice that, while keeping constant all the other param- 
eters , no soliton decay is observed when the disorder 
potential amplitude is lowered to zero. 

Other similar experiments have been performed for 
mesas of different dimensions with varying initial flow 
velocity (see Fig.3S in online supplementary materials). 



In agreement with a previous work [20[ , vortex pairs are 
nucleated for smaller injection momentum. On the other 
hand, an irregular regime with the nucleation of many 
quantized vortices is observed for a mesa with a diameter 
of 9 jim (see Fig.4S in online supplementary materials), 
as predicted for atomic BEC [361 ]. 

In conclusion, we have observed, in a dissipative po- 
lariton gas, oblique dark solitons resulting from the scat- 
tering of a polariton wavepacket against an engineered 
circular attractive obstacle. Soliton properties can be 
controlled by tuning the initial density and velocity of 
the fluid. Stability conditions has been quantitatively in- 
vestigated due to the time resolved polariton phase and 
amplitude. Solitons tend to reach the maximum possible 
value of the ratio n s / no within a time that gets shorter for 
higher excitation density. At such a point, cavity disor- 
der was found to be responsible for the soliton instability 
and their decay into vortex streets. 
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